Effective spontaneous VT -symmetry breaking in hybridized metamaterials 
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We show that metamaterials can be used as a testing ground to investigate spontaneous symmetry breaking 
associated with non-Hermitian quantum systems. By exploring the interplay between near-field dipolar coupling 
and material absorption or gain, we demonstrate various spontaneous breaking processes of the VT -symmetry 
for a series of effective Hamiltonians associated to the scattering matrix. By tuning the coupling parameter, 
coherent perfect absorption, laser action and gain-induced complete reflection (71 reflector) by using an ultra- 
thin metamaterial can be obtained. Moreover, an ideal ?*7~- symmetry can be established effectively in a passive 
system by using metamaterials. 

PACS numbers: 78.67.Pt,42.25.Bs,78.20.-e 



Metamaterials, periodic arrangements of artificial subwave- 
length atoms, have paved the way for realizing different novel 
concepts and applications in optics, such as negative re- 
fraction, sub-diffraction imaging and invisibility cloakingljj- 
0]. Instead of using the macroscopic material parameters 
for description, many of the optical properties, including 
metamaterial-based chirality[5, 6] and electromagnetically in- 
duced transparency ( HIT') lr7 M14f1 . can be understood by study- 
ing the coupling between different magnetic and electric artifi- 
cial atoms in the so called hybridization picture 



Similar to optical lattices with controllable coupling] 17- 
the design flexibility in the coupling between artificial atoms 
opens up a new perspective to realize quantum simulators 
to achieve effects which are difficult to obtain in quantum 
svstems fll7l Il8ll . Here, we theoretically investigate the pos- 
sibility in constructing a metamaterial system with an effec- 
tive VT -symmetric Hamiltonian which can still show real 
eigenenergies. One of the most interesting properties of a 
VT -symmetric Hamilton is the existence of a phase tran- 
sition arising from spontaneous symmetry breaking by tun- 
ing an external parameter, defined as the coalescence of two 
eigenvalues and their corresponding eigenvectors [20]. Un- 
like the previously considered symmetry breaking effects of 
metamaterialsyjj, such a spontaneous symmetry breaking is 
implicit in nature. It means that the phase transition occurs 
without structurally destroying a specific symmetry (e.g. mir- 
ror, rotational symmetry). The external parameter being tuned 
still respect the VT -symmetry. 

During the past few years, VT -symmetric Hamiltoni- 
ans have in fact been extensively investigated within the 
framework of optics. The most direct analog is repre- 
sented by a pair of optical waveguides with a balanced gain 
and loss profile lfl9l H^-fUl- Spontaneous VT -symmetry 
breaking has also been theoretically and experimentally 
investigatedjljjl 22-31]. Intriguing phenomena, such as 
asymmetric transmissi onlfl9[ l24ll . power oscillation ll24ll . uni- 
directional invisibility 4l9l |29[], coherent perfect absorption 
and laser ll3oTl have been predicted. Moreover, extending from 
an ideal !PT-symmetric system with balanced gain and loss, 
non-ideal VT -phase transition of a passive system has also 



been proposed and experimentally demonstrated 125L 12711 . 
Unconventional loss-induced optical transparency is demon- 
strated around the phase transition point[25]. Actually, VT - 
symmetry has not been discussed within the context of meta- 
materials. By employing metamaterials, the larger flexibility 
in designing the coupling between different artificial atoms 
gives us additional possibilities in realizing non-trivial VT - 
symmetry. Here, we show that by mapping a metamaterial 
consisting of a bright electric and a dark magnetic atom to 
an effective Hamiltonian, we can establish an ideal sponta- 
neous symmetry breaking transiting from little absorption to 
complete absorption even in a passive system. It relies on the 
inverse relationship that scattering is much larger than loss for 
a bright electric atom while the reverse is true for a dark atom. 
Such a system has been previously considered to achieve EIT. 
However, we are interested in the phase transition to obtain 
complete absorption Oa 13311 across a threshold in the cou- 
pling parameter, in contrast to the high transmission for EIT 
emerging from the strong coupling regime|7|. Moreover, by 
incorporating a small material gain into the dark atom of the 
ultrathin metamaterial, we can obtain laser action IT34l - [36ll and 
gain-assisted complete reflection by considering the sponta- 
neous symmetry breaking effect of a class of Hamiltonians 
associated with the scattering matrix. 

The discussion of !P7~-symmetric Hamiltonians originates 
from Bender and colleague's work, showing that a Hamilto- 
nian H respecting a parity-time symmetry, i.e. [VT,H] = 0, 
can have completely real eigenenergies [20]. To find the ana- 
log of !PT-symmetry in a metamaterial system, the simplest 
way is to look for a pair of coupled artificial atoms, since 
their response can be understood within a hybridization pic- 
ture in quantum mechanics as a two state problem, and the 
!PT-symmetry in two-state quantum system is already well 
studied[20]. Here, we consider a metamaterial constituting a 
planar array of atoms. Each unit cell consists of an electric 
atom coupled to a magnetic atom through the near fields. For 
simplicity, we assume normal incidence and there is no polar- 
ization conversion for the incoming plane waves. The optical 
response of such a metamaterial can be expressed using the 
scattering matrix S as 
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where a± and b± are the E-field amplitudes of the incoming 
and outgoing waves in the positive and negative directions. 
The same scattering process has a microscopic representation 
using the atomic responses and their associated far-field ra- 
diations. By further assuming each atom is dominated by 
only one dipolar resonating mode, the atomic responses can 
be summarized as 
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while the outgoing fields and the atomic responses are related 
to each other by 
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In the above equations, y ; = y s + y loss {i = 1,2) represents 
the total resonating linewidth for each atom where y s , a pos- 
itive number, represents the part of scattering and y loss is the 
part due to absorption loss when it is positive. We also nor- 
malize the electric dipole p and the magnetic dipole m by the 
oscillator strengths as in Ref.7, and assume the two resonant 
frequencies at the same frequency co hres — coq. In fact, Eq. (4) 
defines a series of Hamiltonian H n with meaning in the fol- 
lowing. By taking an appropriate arithmetic average of Eq.(2) 
and Eq.(3), one can prove that det(H„ - coT) = yields an 
eigenvalue (« - 2)/n of matrix S. We can also obtain 



, det(H 2 - col) det(Hi-cj) 

det(S) = = 2 1 

det(Ho-wI) det(Ho-wI) 
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For example, det(Ho - col) = leads to the pole of S, which 
corresponds to laser action ll30l 34-361. det(H2 - col) = 
implies a zero of S to completely absorb the incoming EM 
waves, namely coherent perfect absorption 
For a metamaterial with mirror symmetry in the propagation 
direction, the two eigenvectors of S are (a+,a_) = (1, 1) with 
eigenvalue t+r and (1, -1) with eigenvalue t—r by defining t/r 
as the complex transmission/reflection coefficient for electric 
field with respect to the mirror plane. If we further confine our 
discussion to the coupling between a bright electric and a dark 
magnetic atom, we have the additional requirement t = 1 + r. 
Then, det(Hi - col) = implies a complete reflection r = — 1 
and t = 0. In general, a particular n/eigenvalue of S denotes 
a particular optical phenomenon and it corresponds to solv- 
ing the eigenvalue problem of H„. If the eigenvalue of H„ 
is real, the phenomenon can occur at a real frequency with a 
pole/zero on the corresponding optical spectrum. Therefore, 



by considering the spontaneous symmetry breaking process 
of a VT -symmetric H„, one can observe the phase transi- 
tion when the real eigenfrequencies are turned into complex 
conjugates or vice versa. For our 2x2 Hamiltonian, the P- 
operation switches rows and columns of the matrix while the 
T-operation takes complex conjugate on the elements. A VT - 
symmetric H n is defined by \PT, H„] = 0. We have to empha- 
size that PT -symmetry being discussed here is an effective 
one through the establishment of an effective !PT-symmetric 
Hamiltonian, which is different to the approach employing 
complex potentials j20ll . 

As an example, a VT -symmetric H 2 corresponds to the 
condition j2 - 2y s 2 = 2y\ - y\, which can be satisfied by a 
representative hybridized metamaterial with a bright electric 
atom (scatters most) and a dark magnetic atom(absorbs most). 
For a fixed set of parameters, e.g. j\ = 15 THz, y'° ss = 1.75 
THz (yj = 13.25 THz), 72 ~ y^ ss , this condition is repre- 
sented by the vertical red line in Fig. 1 at 72 = H-5 THz 
to indicate a VT -symmetric H 2 . Then, the spontaneous sym- 
metry breaking process can be observed by changing the cou- 
pling coefficient k (still respecting the !PT-symmetry) across a 
threshold value. The eigenfrequencies of VT -symmetric H 2 
are real(solid red line) when k > yi, and becomes a com- 
plex conjugate pair(dotted red line) when k < 72- The split- 
ting of the real/imaginary part of the eigenfrequencies with 
varying k is showed in the inset. This behavior reveals itself 
through the absorption coefficient when the metamaterial is 
excited by the incident waves of (a + , a_) = (1,1), giving CPA 
for k > 72 at the two real frequencies. We note that such 
a configuration of coupling between a bright and a dark atom 
(without the ^"T-symmetric condition satisfied in general) has 
been used to give EIT in the strong coupling regime[7J. We 
endow the same system to visualize effective !P7~-symmetry 
breaking and focus on its optical properties in the regime near 
the critical point. In the following, we demonstrate the inter- 
esting optical properties associated with VT -symmetry break- 
ing with metamaterial structures. In particular, undergoing the 
phase transition, the metamaterial becomes perfectly absorp- 
tive with real frequencies. The bright atom (inset of Fig. 2(a)) 
is a silver strip of length LI = 150 nm, width Wl = 50 nm 
and thickness t — 20 nm. Fig. 2 (a) shows the electric field 
near one end of the strip with numerical simulation (symbols) 
and theoretical model (lines). The numerical simulation is car- 
ried out using finite element analysis (Comsol MultiPhysics) 
with Drude model of silver being e = e M - to 2 p {<x) 2 + iy p cS) , 
where to p = 2.1961 x 10 15 Hz, y p = 4.3439 x 10 12 Hz, and 
€ m = 3.7. The theoretical result is fitted using a Lorentz model 
E s cc (coo -co- iyxT\y\ - y[ oss )Eo, where co = 402.5 THz, 
71 = 15 THz, and y'° ss = 1.75 THz, E m is the scatter- 
ing(incident) E field. y ° ss can be estimated by the probing the 
response of the bright atom with and without metal loss. The 
dark atom is a magnetic resonator tuned to the same resonat- 
ing frequency with two parallel silver strips and a dielectric 
gap material (e = 2.25 + ie" with e" introduced to represent a 
slight absorption or gain) in between. The geometric param- 
eters are L2 = 100 nm, W2 = 30 nm, t = 20 nm, as showed 
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FIG. 1 : (Color online) Parametric evolution of the metamaterial for 
different H„(n = 0, 1,2); Inset: Splitting of eigenvalues for H2 with 
varying k. Vertical solid lines indicate the real eigenvalues case for 
ideal VT -symmetric H n above the phase transition threshold, short 
dotted lines illustrate the complex conjugate pair eigenvalues case for 
ideal VT -symmetric H n below the phase transition threshold. Inset 
gives the eigenvalues with varying k for VT -symmetric H2. Curved 
solid lines show the real eigenvalue condition indicating the optical 
properties for non-ideal VT -symmetric H„, short dashed lines show 
branching condition for non-ideal VT -symmetric H n . 



in the inset of Fig. 2(b). Due to the large quality factor of 
the dark mode resonance, a small e" of the gap material can 
effectively alter the half-linewidth 72. As shown in Fig. 2 (b), 
numerically calculated field (symbols) near one of the metal 
strip ends is again fitted well with theoretical Lorentz model 
(solid line) for e" = 0,-0.1 (hollow, solid symbols). Com- 
paring between numerical and analytical results with different 
e", a simple approximated formula between 72 and e" can be 
established as 72 = 1-75 + 22. 5e". Now, we hybridize these 
two atoms within a metamaterial with a square lattice of pe- 
riodicity d = 400 nm and a distance h in separating the two 
atoms, as shown in Fig. 2(c). The change of distance h be- 
tween the two meta-atoms can tune the coupling k between 
them, but there exists no simple relationship between them, as 
shown in Fig. 2(d). The k can be calculated from the bright 
meta-atom probe response for e" = by changing the dis- 
tance h between them. The fitted formula can represent main 
variation character of k, but exhibits slight deviation in large k 
region. 

Now, we illustrate the spontaneous symmetry breaking with 
fixed geometry of the meta-atoms with constant a>o, 71, and 
yloss fjjg variable parameters are the coupling coefficient k 
(or equivalently h) between the two atoms and 72 (through 
the tuning of e"). We can then choose a particular 72 to sat- 
isfy !PT-symmetry requirement of H„. For a VT -symmetric 
H 2 , we need 72 = 27; - 7, = 11.5 THz, i.e. e" = 0.433, 
which is a passive configuration to display the effective VT 
symmetry. In the light of the above theoretical discussion, the 
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FIG. 2: (Color online) Analytical (solid line) and numerical (scatter 
line) spectral response of the E-field probe for bright(a) and dark(b) 
meta-atom to establish the atomic Lorentz model with fitted param- 
eters (a> THz, y! THz, y'° ss , and y 2 ). (c) and (d): The hybridized 
metamaterial under consideration and the relationship between the 
coupling coefficient k and the distance h between the two meta- 
atoms. 
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FIG 3: (Color online) Observation of effective spontaneous VT- 
symmetry breaking for H 2 and Hi. Analytical model (with fitted 
parameters from Fig. 2) (a) and numerical(b) outgoing \t + r\ 2 spectra 
versus frequency and coupling distance. Analytical(c) and numeri- 
cal(d) transmission spectra versus frequency and coupling distance. 

eigenvalue problem for a !P7~-symmetric H2 is equivalent to 
the zeros of the S, which can be easily observed in the co- 
herent outgoing \r + t\ 2 spectra when we change the param- 
eter h. The VT -phase transition point in this case occurs at 
Kpj = 11.5 THz (h = 46.0 nm). Below the critical point, 
Kpr < 11.5 THz (h > 46.0 nm), the eigenvalues of H2 is 
a complex conjugate pair, there are no crossing in real fre- 
quency axis, so we cannot observe CPA in this region, shown 
in Fig. 3(a). As we increase (decrease) k(Ii), at Kpt = 11 -5 
THz(/i = 46.0 nm), the two eigenvalues of H2 collapse to 
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one point in the real frequency axis, i.e. CPA(|r + t\ — 0) 
emerges in the spectra. When k passes over the Kpj, there 
are two real frequencies satisfies the eigenvalue solution of 
H2, CPA branching phenomenon can be clearly shown in Fig. 
3(a) from the analytical results obtained from Eq. (2) and (3) 
with the atomic model established from the previous simula- 
tions shown in Fig. 2. The same phenomenon can be ob- 
served from the numerical simulations (with microstructures 
included) shown in Fig. 3(b). Real parts (white dashed line in 
Fig. 3(a), (b)) of the eigenvalues of H2 agree well with the CPA 
positions in the spectra. In both Fig. 3(a) and Fig. 3(b), the 
absorption loss (1 - \t + r\ 2 for symmetric incoming waves) of 
the hybridized system increases dramatically with increasing 
coupling (decreasing separation) below Kpj. At Kpj, the pres- 
ence of m interferes with p to the right size and gives r = — I 
and t = A. This is the onset of CPA phase. Going beyond 
the transition point (larger coupling), the CPA peak splits into 
two. On the other hand, this splitting makes both atoms can- 
not be efficiently excited at a>o and loss decreases at that fre- 
quency. In fact, the non-monotonic function of loss against 
coupling coefficient is a direct manifestation of spontaneous 
symmetry breaking, which is similar to the loss-induced opti- 
cal transparency in a passive !PT-symmetric system deviated 
from the ideal VT -symmetry B25L Here, using metamaterials, 
we can actually establish an ideal !PT-symmetry in a com- 
pletely passive system. CPA becomes possible besides the 
counterintuitive trend of loss function. 

Besides H2 in revealing !PT-symmetry breaking process, 
Hi is also a potential candidate to display the same process. 
As predicted in Eq. (5) and Fig. 1, by taking into account 
the eigenvalue problem of Hi, we need to fulfill the case 
72 = -y'° ss = -1.75 THz (the black vertical line in Fig. 1), 
i.e. e" = -0.156, where potential amplifying media can be 
available in nature 1 34- 3^1 . The emergence of effective VT- 
symmetry breaking process can be observed though transmis- 
sion spectra, and t = 0(r = -1) in spectra implies real fre- 
quency solution for eigenvalue of Hi. It indicates the system 
without absorption loss, A — 1 - \t + r\ 2 = 0, and transmis- 
sion, t — 0, with perfect reflection with and a 7r-phase shift, 
r = — 1. The phase transition point in this configuration is at 
Kpt = 1.75 THz (h = 98.4 nm). Illustrated in Fig. 3(c), and 
numerically verified in Fig. 3(d), the phase transition point 
is at h = 97.5nm, t = (r = -l)in the transmission spec- 
tra agrees well with the real parts of the eigenvalues of Hi. 
Below the phase transition point (k < »pr), the eigenvalue of 
H\ implies no cross point in the real frequency axis, a com- 
plex conjugate pair, so transmission is dominated by the bright 
atom dipole oscillation, although f * 0, the eigenvalue of S is 
not -1 (Absorption loss is not zero). Above the phase transi- 
tion point Kpj-, the eigenvalue of H\ implies two real eigen- 
frequencies, so t = (r = — 1) can be observed at two differ- 
ent frequencies in the transmission spectra. Real parts (white 
dashed line in Fig. 3(c), (d)) of the eigenvalues of Hi match 
the t = positions well in spectra. The effective spontaneous 
VT -symmetry breaking for Hi can thus be observed in the 
spectra. It is worth to note that the transmission exhibits maxi- 



mum value at the position k s = y-7271 = 5.12 THz (h = 73.8 
nm for analytical prediction and h = 72.5 nm for numerical 
results) at ojq, because gain behavior of the system appears, 
when pass over the phase transition point Kpr- As we in- 
crease k, we can encounter a real eigenvalue of Hq, indicating 
pole of S, which implies t he sy stem shows laser action at this 



parameter position! 30, 34-36], as indicted in the cross point 



between the black vertical line indicating !PT-symmetric Hi 
and the blue curve displaying real eigenvalue of Hq in Fig. 1. 
The laser action results from non-ideal VT -symmetry of Hq 
whose origin will be discussed in the following. 

For Ho, ideal !PT-symmetry condition 72 = -71 implies 
that huge gain is required to compensate metal losses. This 
is difficult to be satisfied using natural materials. However, 
in a practical (non-ideal) system, even though ideal condi- 
tion is not satisfied, the eigenvalues still split in the same way 
as an ideal !PT-symmetric Hamiltonian. Only an additional 
constant bias in the imaginary part of the eigenvalues occurs. 
Similar strategies have been employed in studying non-indeal 
VT -symmetry in passive svstems ll25l, 12711 . Then, phase tran- 
sition (branching) point is not on the real frequency axis, and 
is shifted into the lower-half of the complex w-plane. One of 
the eigenvalues (the one with increasing imaginary part) will 
cross the real frequency axis, corresponding to the original 
optical phenomenon, e.g. laser action for the case of Ho as a 
pole of S at k s - y- yi yp IT30l 34 - 3^] . It is the solid blue curve 
in Fig. 1 and is lying below branching condition (branching 
occurs when k > Kpj-) where k<pt = (71 - y2)/2(dashed green 
line). In contrary, ideal TT -symmetric Ho has the poles of S 
on the real frequency axis beyond branching condition (ver- 
tical solid blue line) at two different frequencies. As an ex- 
ample, 72 = -0.5 THz (e" = -0.1), Fig. 4 gives the outgo- 
ing intensity of the system with symmetric incoming light as 
a function of k(Ii). Laser action collapses to one position at 
u>o, h = 89.4 nm for theoretical prediction and h = 90.0 nm 
for numerical results. The branching phenomenon still occurs 
(without real poles of S) and agrees well with the real parts 
of the eigenvalues of Ho (black dashed lines), as displayed 
in Fig. 4. Laser action and branch phenomenon is the direct 
demonstration of non-ideal !PT-symmetry for Ho. Besides 
Ho, non-ideal VT -symmetry can also occur for Hi and H2. 
In Fig. 1, the real eigenvalue of H2 (zero of S), i.e. CPA, 



is at k 



27', '") (solid red curved line) while the 



branching condition is at - (71 + 72 - 7j OS 0/2 (dashed 
purple line). For Hi, perfect reflection, r = — 1, t = (one S 

eigenvalue is —1) is at /c_i = y-727^" (solid black curved 

line) while the branching condition is at K<pj- = {y^ ss - yi)l1 
(dashed black line). So, for a non-ideal !P7~-symmetric sys- 
tem, the optical phenomenon (e.g. laser, CPA or ^--reflector) 
associated to H n collapses to a single point (with a k a bit less 
than the one triggers branching) instead of the two branches. 

In the above discussion, k is the external parameter in 
demonstrating !P7~-symmetry. If we change 72 with a fixed 
k instead (travelling in a horizontal direction in Fig.l), the 
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FIG. 4: (Color online) Output intensity (color map) of the metama- 
terial log w \t + r\ 2 revealing spontaneous symmetry breaking of non- 
ideal VT -symmetry for Ho and the real part of the eigenfrequencies 
for Ho(dashed lines), (a) Analytical model; (b)Numerical simula- 
tions with metamaterial microstructures. 
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FIG. 5: (Color online) Analytical (solid line) and numerical (dashed 
line) outgoing spectra, as a function of y 2 f° r the metamaterial at a> . 



metamaterial can function as a laser, the cross point for Ho, 
or a CPA, the cross point for H2 at ojq. As an example for 
demonstration, h — 75 nm, besides appearance of laser ac- 
tion and CPA in spectra, schematically shown in Fig. 5, above 
CPA position, an increase of loss (increasing jj) leads to more 
output of the system, which is similar to loss induced optical 
transparency in a non-ideal passive VT svstem ll25ll . Similar 
counter-intuitive trend also occurs beyond the laser position, 
an increase in gain (decreasing negative 72) causes less output 
observed in the spectra. 

In conclusion, we investigated the effective !PT-symmetric 
Hamiltonians using a metamaterial system to reveal sponta- 
neous symmetry breaking processes. Theoretical and numer- 
ical results clearly confirm the analog of spontaneous break- 
ing of VT -symmetry using metamaterials. Our results not 
only open up a new perspective towards testing the intrigu- 
ing properties of spontaneous !P7~-symmetry breaking in non- 
Hermitian quantum systems, but also provide feasibility for 
future applications in plasmonic and metamaterial systems, 
such as laser actionll30lf34 - 36 ] , and CPA with sub-wavelength 
thickness 032ll . Moreover, by using metamaterials, it becomes 
possible to establish an ideal VT -symmetry by using a passive 
system. 

This work was supported by Hong Kong Research Grants 
Council (GRF grant CityU 102012). 
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